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ABSTRACT 


The  modulated  beam  of  the  free  electron  laser  when 
passed  through  a  thin  target  should  produce  strong 
transition  radiation.  The  measurement  of  the  transition 
radiation  is  a  direct  measurement  of  the  electron  beam 
modulation  in  the  free  electron  laser.  A  transition 
radiation  experiment  using  the  Stanford  MKIII  Infrared  Free 
Electron  Laser  (IRFEL)  has  been  proposed.  The  analysis  has 
centered  on  TRANSPORT,  a  computer  program  used  for  designing 
charged  particle  beam  systems.  The  MKIII  IRFEL  wiggler  exit 
bending  magnet  system  was  modeled  using  TRANSPORT.  Analysis 
reveals  that  the  transverse  emittance  and  momentum  spread 
characteristics  will  cause  the  modulated  beam  to  demodulate 
along  the  path  of  the  central  trajectory.  A  detector 
location  10  millimeters  downstream  of  the  first  bending 
magnet  is  found  suitable  for  the  measurement  of  the  electron 
beam  modulation.  For  this  case  the  thin  foil  must  be 
rotated  approximately  11.25  degrees  about  the  y-axis  to 
yield  a  minimum  effective  picobunch  extent  as  seen  by  the 
thin  foil.  Alternate  FEL  systems  are  suggested  for 
investigation . 
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I.   INTRODUCTION 

A.   OVERVIEW 

The  electron  beam  from  a  free  electron  laser  (FEL) 
offers  the  possibility  of  observing  and  using  transition 
radiation  (TR)  at  optical  frequencies.  After  interacting 
with  the  wiggler  and  radiation  fields,  the  electron  beam 
becomes  modulated  in  the  axial  direction  at  the  lasing 
frequency  [Ref.  1]  .  The  axially  modulated  bunches 
(picobunches)  could  be  passed  through  a  thin  foil  producing 
the  surface  effect  TR  as  illustrated  in  Figure  1-1. 
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Figure  1-1   Transition  Radiation  (TR) 


The  axially  modulated  beam  picobunches  impact  a  thin  foil  of 
aluminum  or  mylar  yielding  forward  and  backward/ specular  TR 
[Ref.  2].  Both  forms  of  radiation  should  be  enhanced  by 
coherent  emission  from  all  charges  in  the  picobunch.  The 
backward/specular  TR  is  redirected  to  a  viewing  port/optical 
detector  by  a  mirror.  The  resulting  radiation  could  be  used 
to  diagnose  the  electron  beam  modulation  [Ref.  1] . 

The  TR  power  radiated  per  picobunch  is  found  to  be 
proportional  to  the  charge  of  the  picobunch  sguared.  For  a 
picobunch  of  approximately  106  electrons,  the  radiation 
intensity  can  easily  be  detected  and  measured  with  an 
optical  detector.  The  axially  modulated  electron  beam  of 
the  IRFEL,  upon  impact  with  the  thin  foil  or  multifoil,  will 
generate  coherent  optical  transition  radiation.   [Ref.  2] 

B.   MODULATION  OF  THE  ELECTRON  BEAM 

The  longitudinal  bunching  of  the  electron  beam  by  the 
undulator  of  the  free  electron  laser  (FEL)  is  a  result  of 
the  force  acting  on  the  individual  particles  by  the  electric 
field  of  the  superimposed  laser  radiation.  This  effect  on  a 
single  electron  is  examined  in  Figure  1-2  [Ref.  3].  This 
figure  illustrates  the  fields  and  trajectories  within  the 
undulator  of  a  FEL:  an  alternating  magnetic  field,  a  co- 
propagating  laser  electric  field  and  an  electron  along  a 
trajectory.  The  electric  field  of  the  laser  is  zero  at 
position  A.  The  electron  experiences  no  force  other  than 
that  of  the  undulator 's  magnetic  field,  which  will   bend   it 
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Figure  1-2   Longitudinal  Bunching  of  the  Electron  Beam 


in  the  horizontal  plane.  Now  consider  the  electron  a  time 
later  at  position  B.  The  electron  has  traveled  one-quarter 
of  an  undulator  period  and  now  experiences  a  maximum  in  the 
laser  electric  field.  The  force  on  the  electron  at  point  B 
is  opposite  the  direction  of  the  laser  electric  field  at 
that  point  and  is  labeled  Fe.  The  electron  will  therefore 
experience  a  decelerating  force.  Whether  an  electron  is 
accelerated  or  decelerated  by  the  laser  electric  field 
depends  on  the  phase  of  the  laser  electric  field  with 
respect  to  the  electron  oscillation.  Some  electrons  are 
accelerated,  while  those  half  a  laser  wavelength  ahead  or 
behind  are  decelerated.  For  example,  an  electron  at  C  would 
be  accelerated.   As  the  accelerated  electrons  catch  up  to 


those  which  have  been  decelerated,  they  form  into  bunches  at 
the  laser  wavelength.  In  other  words,  the  electron  beam 
becomes  axially  modulated  at  the  laser  wavelength  of  the  FEL 
radiation.   [Ref.  3] 

Figure  1-3  [Ref.  4]  depicts  approximately  106  electrons 
being  bunched  within  the  optical  wavelength  of  the  FEL 
radiation  while  within  the  undulator. 
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Figure  1-3   Bunching  of  Electrons  within  the  Wavelength 
of  Laser  Radiation 


As  the  electrons  oscillate  they  spontaneously  emit  "magnetic 
bremsstrahlung"  peaked  in  the  forward  direction  thus 
generating  the  laser  electric  field.  When  the  laser 
frequency  is  nearly  resonant  with  the  electron  oscillations, 
significant  amounts  of  energy  can  be  exchanged  between  the 
electrons  and  the  laser  field  causing  the  bunching  effect  of 
the  electron  beam.  The  condition  for  this  resonance  is 
given  by  Equation  1-1  (MKS)  [Ref.  5]: 


y2    -    AW/2A(1   +   e2B2ms  V2/4 TT2m02c2 )  ,  (1-1) 

where : 

Y  =  energy  of  resonant  electron  divided  by  rest 
mass  energy, 

\j  =  wiggler  (undulator)  magnetic  period, 

A  =  laser  wavelength  (operating  wavelength) , 

e  =  electron  charge, 

Brms  =  rms  magnetic  induction  of  the  wiggler, 

c  =  velocity  of  light  in  vacuo. 

The  operating  (lasing)  frequency  of  a  FEL  is  determined 
by  Equation  1-2  (MKS) , 

A  =  Aw/2y2  (1  +  K2) ,  (1-2) 

where : 


K  =  the  dimensionless  measure  of  the  wiggler 
strength, 


K  =   eBrmsAw/27Tmoc-  C1"3) 

The  operating  parameters  of  the  Stanford  MKIII  IRFEL  are 
shown  in  Table  1-1  [Ref.  6]. 

The  wiggler  period  of  the  IRFEL  is  of  the  "hybrid"  type 
employing  both  SmCo  permanent  magnets  and  vanadium  permendur 
steel  poles.   The  magnets  are  oriented  to  provide   a   nearly 


TABLE  1-1 
MKIII  IRFEL  PARAMETERS 


WIGGLER 
magnet  period 
number  of  periods 
minimum  gap 
wiggler  strength 
peak  magnetic  field 
RESONATOR 


N 

9min 

K 

B^ 


2.3 

47.5 

6.5 

1.14 

0.702 


cm 


mm 


length 

L 

178.57 

cm 

Rayleigh  range 

ZR 

73 

cm 

E-BEAM 

energy 

Y 

85 

energy 

E 

43 

MeV 

energy  spread 

<5yA 

0.5 

% 

peak  current 

I 

20 

Amps 

emittance 

£N 

7 

it— mm  mrad 

macropulse  length 

T 

3 

ysec 

sinusoidal  vertical  field  By,  using  a  Halbach  undulator 
design.  [Ref.  6]  The  lasing  frequency  can  be  modified  by 
slightly  changing  the  energy  of  the  beam  or  by  changing  the 
K  value  of  the  wiggler.  In  practice  the  K  parameter  is 
modified  by  adjusting  the  gap  between  the  permanent  magnets 
in  the  wiggler  cavity.  The  resonator  length  must  then  be 
modified  to  account  for  the  new  wavelength.   Assuming  the 
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lasing  wavelength  of  3.0  microns,  the  IRFEL  will  produce  a 
modulated  beam  as  shown  in  Figure  1-4. 
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Figure  1-4  Disks  of  Charge  (picobunches) 

The  rectangles  indicate  disks  of  charge  of  thickness  1 
micron  separated  by  3  microns.  This  periodic  bunched  pulse 
train  will  generate  transition  radiation  of  wavelength  3 
microns  and  also  odd  harmonics  when  passing  through  a  thin 
foil  such  as  mylar  or  aluminum.  The  current  vs  time  profile 
is  shown  in  Figure  1-5. 


Figure  1-5   Current  pulses 


C.   WHAT  IS  TRANSITION  RADIATION 

Transition  radiation  (TR)  is  the  electromagnetic  radiation 
which  is  emitted  when  a  uniformly  moving  charged  particle 
traverses  from  one  medium  into  another  with  a  different 
dielectric  constant.   More  generally  the  effect  will  occur 
in  the  presence  of  a  heterogeneity  in  a  medium.   [Ref.7] 

TR  may  be  considered  to  originate  from  the  change  in  the 
fields  as  an  electron  crosses  a  boundary.  The  magnetic 
field  remains  essentially  constant  but  the  electric  field 
changes  because  of  the  change  in  the  dielectric  constant.  A 
third  field,  radiation,  is  required  to  meet  the  boundary 
conditions. 

TR  can  also  be  explained  in  terms  of  a  charge  and  its 
image.  While  an  electron  moves  uniformly  in  vacuo  it  does 
not  radiate;  the  field  in  vacuo  is  equal  to  the  field  of  the 
electron  and  its  image  moving  toward  it.  When  the  electron 
crosses  the  boundary  surface  both  the  electron  and  its  image 
cease  to  exist  from  the  point  of  view  of  the  field.  The 
crossing  of  the  boundary  must  give  rise  to  exactly  the  same 
radiation  as  that  which  is  due  to  a  sudden  stop  at  the  same 
point  of  the  electron  and  its  image.   [Ref.  7] 

The  free  electron  laser  will  longitudinally  modulate  the 
electron  beam  in  the  undulator/wiggler  yielding  disks  of 
charge  of  small  thickness  (lu) .  Due  to  the  low  efficiency 
of  OTR,  106  electrons  per  picobunch  will  yield  approximately 
104  photons  per  picobunch  upon  foil  impact. 

For  a  relativistic  beam  of  electrons  (  g  ->  1)  ,  and  a  TR 
foil  inclined  at  45  degrees  with  respect  to  the  beam  axis 


(Figure  1-1)  ,  the  backward/specular  reflection  intensity  for 
a  vacuo/dielectric  interface  is  given  by  the  relation  [Ref. 
2]: 

d2I/dwdft  =   F(  i|j,w)  (e2/^2c)  [62/(y~2  +  82)2]  ,       (1-4) 

where : 

,     2/w.lQN1l/2    _    „    _„/fIlJ.ft,     2 


„ ,  ,  r  [e-sin*(i|H-9)]     7        -    £    COSjjp+ejU 

F(^/w)  =  [ 2 172 TJ 

[e-sin  (ifH-6)J  7   +  e  cosOp+6) 


F  =  Fresnel  coefficient  for  reflection  of  light 
with  the  E-vector  in  the  plane  of  incidence, 

ijj  =  angle  of  specular  reflection  (45  degrees)  , 

£    =   dielectric  constant  of  foil  material. 


The  intensity  of  this  radiation  should  be  sufficient  for 
the  detection  and  observance  of  TR  from  a  single  foil/multi- 
electron  interface.  Measurement  of  the  radiation  will  be 
accomplished  via  a  detector  of  the  phototube,  TV  camera,  or 
reticon  type. 

The  characteristics  of  TR  that  will  apply  to  the 
proposed  experiment  are: 

(1)  The  angular  distribution  of  TR  is  such  that  for 
relativistic  particles  the  emission  angle  of  the 
photons  with  respect  to  the  electron  beam  is 
appproximately  equal  to  1/y . 

(2)  The  intensity  of  TR  is  a  strong  function  of  the 
energy  of  the  particle  producing  it. 


(3)  The  radiation  is  polarized  with  the  electric  field  in 
the  plane  containing  the  particle  velocity  vector  and 
the  wave  propagation  vector. 

(4)  In  the  case  of  a  periodic  medium  (i.e.,  multifoil)  an 
rtr  phase  difference  is  the  condition  for  constructive 
interference  or  resonance  between  radiation  produced 
at  successive  interfaces;  r  is  a  positive  integer. 
For  m  boundaries  the  intensity  varies  as  m2 . 

(5)  If  the  medium  is  a  transparent  dielectric  then 
e'=n2(w),  where  n  is  the  refractive  index  and  e1  is 
the  dielectric  constant. 

(6)  For  a  relativistic  electron,  TR  will  arise  at  the 
surface  of  the  interface.   [Ref.  7] 

(7)  The  power  radiated  from  a  transition  foil  is 
proportional  to  the  charge  squared  [Ref.  2]. 
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II.   PROPOSED  EXPERIMENT  USING  MKIII  IRFEL 

A.   EXPERIMENT  PROPOSAL 

The  purpose  of  this  experiment  will  be  to  measure  the 
electron  beam  modulation  of  the  FEL  using  a  TR  foil. 

The  MKIII  IRFEL,  located  at  the  High  Energy  Physics  Lab 
at  Stanford  University  may  lend  itself  to  an  experiment  to 
measure  the  axial  modulation  of  the  electron  beam  of  the 
FEL.  The  "bunched"  beam  downstream  of  the  wiggler  magnet  is 
normally  bent  away  from  the  downstream  cavity  mirror  and 
directed  to  the  beam  dump.  If  the  beam  were  to  retain  its 
axial  modulation,  upon  the  passage  through  a  thin  foil  it 
would  generate  coherent  transition  radiation.  The  thin  foil 
will  hereafter  be  referred  to  as  the  detector.  The  detector 
must  be  conveniently  located  downstream  of  the  wiggler. 

1.   System  Description 

The  MKIII  IRFEL  is  shown  in  Figures  2-1,  2-2,  2-3, 
and  2-4.  Figure  2-1  is  a  photo  of  the  accelerator  section 
of  the  IRFEL.  Figure  2-2  is  a  representation  of  the  IRFEL 
downstream  of  the  accelerator.  The  electron  beam  produced 
by  the  accelerator  is  directed  around  the  upstream  optical 
cavity  mirror  and  through  the  wiggler  magnet.  The  beam  is 
axially  modulated  within  the  wiggler  and  is  then  directed 
around  the  downstream  movable  mirror  to   the   beam  dump.   An 
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Figure  2-1   MKIII  IRFEL  Accelerator 
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external  photo  of  the  MKIII  IRFEL  wiggler  magnet  is  shown  in 
Figure  2-3. 


Figure  2-3   MKIII  IRFEL  Wiggler 

Figure  2-4  is  a  photo  of  the  experimental  beamline  apparatus 
downstream  of  the  wiggler  cavity.  The  electron  beam  will 
proceed  from  left  to  right.  The  figure  shows  the  waveguide 
type  piping,  flanges,  bending  magnet  #2,  bending  magnet  #3, 
and  a  four  way  cross.  Bending  magnet  #1  is  upstream  to  the 
left  of  the  flange  and  bellows  (not  shown  in  the  photo)  . 
The  detector  could  easily  be  placed  on  the  waveguide  piping 
or  at  the  four  way  cross  location. 
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Figure  2-4   MKIII  IRFEL  Wiggler  Exit 

2 .  Possible  Detector  Locations 

Figure  2-5  is  a  top  view  representation  depicting 
possible  locations  for  the  detector.  The  electron  bunches 
exiting  from  the  wiggler  are  deflected  by  the  first  bending 
magnet  in  order  to  prevent  them  from  striking  the  resonator 
mirror. 

3 .  Hardware  Design  at  Detector  Locations 

The  hardware  design  needed  for  a  detector  between 
locations  B  and  C  is  shown  in  Figure  2-6.  Placement  of  the 
detector  between  points  B  and  C  would  require  modifying  the 
waveguide  piping,  bellows,  and  flange.  Location  E  is  the 
alternative  location  of  the  detector  and  would  require  the 
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Figure  2-5   Detector  Locations 
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Figure   2-6      Detector   Placement   on  Waveguide 
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modification  of  Figure  2-7A.  The  linear  motion  feedthrough 
needs  to  be  connected  to  a  sliding  sleeve  device  in  order  to 
position  the  TR  foil  and  mirror.  A  suitable  sliding  sleeve 
device  is  shown  in  Figure  2-7B.  In  general,  it  would  be 
easier  to  modify  the  IRFEL  by  placing  the  detector  at 
location  E. 

B.   ANALYSIS  OF  ELECTRON  BEAM  TRANSPORT  EQUATIONS 

Note:  Readers  not  familiar  with  the  properties  of  beams  and 
the  TRANSPORT  development  of  system  equations  should  read 
Appendices  A  and  B  before  proceeding  with  this  section. 

The  development  of  the  electron  beam  transport  equations 
for  locations  A  through  E  along  the  central  trajectory  is 
found  in  Appendix  B.  It  is  desired  to  find  the  best 
location  in  the  beamline  to  place  the  detector  in  order  to 
generate  coherent  transition  radiation.  Of  particular 
concern  is  the  longitudinal  extent  of  the  picobunches 
(Figure  1-3)  as  they  traverse  the  beamline.  In  the  ideal 
case  the  picobunches  would  retain  their  axial  modulation 
traversing  the  beamline.  This  is  the  case  of  a  non- 
divergent,  monoenergetic  beam.  In  the  real  world,  electron 
beams  are  non-monoenergetic  and  are  divergent.  Beams  of 
this  nature  are  specified  in  terms  of  their  phase  space.  In 
order  to  observe  coherent  transition  radiation,  the 
picobunches  must  not  "debunch"  significantly  upon  traversing 
the  beamline.  "Significant"  debunching  for  this  experiment 
is   on   the   order   of   a   few   (>3)  microns.   If  debunching 
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exceeds  3  microns  the  coherent  TR  effect  will  not  be 
observed.  A  simple  representation  of  the  beam  traversing 
the  system  is  desired.  Figure  2-8  is  a  representation  of 
the  paths  a  single  electron  located  at  "a"  may  take  in 
traversing  the  system.  System  parameters  (obtained  from 
blueprint)  and  definitions  applicable  to  the  representation 
are  shown  below  the  figure. 


*-  bending  magnets 
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L  centi 


*  x 
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Figure  2-8   Trajectory  of  a  Single  Electron 


x0  =  the   initial   horizontal   displacement  of  the 
electron  with  respect  to  the  central  trajectory. 

0O  =  the  divergence  of  the  electron  with  respect  to 
the  assumed  central  trajectory. 

60   =  the   fractional   momentum   deviation   of   the 
electron  from  the  assumed  central  trajectory. 
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The  central  trajectory  parameters  are: 

Lx  =  drift  length  1  327.03  mm 

L2  =  drift  length  2  9  5.30  mm 

L  =  effective  pathlength  through  3  8.35  mm 

bending  magnet 

a  =  bending  angle  11.25  ° 

Po  =  bending  radius  of  magnet  195.30  mm 

L0  =  initial  drift  length  0 . 0  mm 

E  =  energy  of  electron  40  MeV 

(Note:  L0  =  0;  the  electron  starts  at  the  magnet  pole  face) 

1.   Buskirk  System  Equations 

The  system  equations  for  this  model  were  derived 
by  Dr  Fred  Buskirk.  These  equations  assume  the 
initial  drift  length,  L0  and  the  effective 
pathlength  through  the  bending  magnet,  L,  equal 
zero.   The  derived  equations  are: 

dx  =  (L1/cos2a  +  L2)6  0  +  (L1sina/cos2a)  6  Q  +  xQ,    (2-1) 

ds  =  -(L1sina/cos2a)9 0  -  (L1tan2a)6Q,  (2-2) 


where : 


dx  =  the  change  in  the  horizontal  distance  of  the 
electron  with  respect  to  the  central  trajectory, 
and 

ds  =  the  difference  in  the  pathlength  between  the 
arbitrary  ray  and  the  central  trajectory. 
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The  electron  is  to  traverse  the  system  from  point  "a"  to  the 
xy  plane  of  location  E.  Assuming  the  fringing  fields  are 
negligable  the  electron  will  bend  in  the  xz  plane  only. 

Trajectory  1  (Tl)  is  for  the  case  of  the 
monoenergetic  ( 50  =  0)  non-divergent  (SQ  =  0)  electron 
initially  located  at  x0.  Its  trajectory  coincides  with  the 
central  trajectory  as  expected.  Trajectory  2  (T2)  is  the 
case  of  a  non-monoenergetic  non-divergent  electron  initially 
located  at  x0.  Since  this  electron  possesses  greater 
momentum  than  the  central  momentum  it  will  not  be  bent  as 
much  by  the  bending  magnets.  The  pathlength  between  the 
bending  magnets  will  be  shorter  than  L^.  If  two  electrons 
having  equal  velocities  were  located  at  "a"  start  downstream 
on  trajectories  Tl  and  T2 ,  the  electron  on  T2  would  have  a 
shorter  pathlength  through  the  system.  At  the  time  when  the 
electron  on  Tl  reached  the  xz  plane  of  location  E,  the 
electron  on  T2  would  have  already  passed  throught  the  system 
by  a  distance  equal  to  the  pathlength  difference.  T2  is 
also  slightly  greater  than  Tl  in  the  x  direction. 
Trajectory  3  (T3)  is  the  case  of  a  non-monoenergetic, 
divergent  electron  initially  located  at  x0.  As  this 
electron  traverses  the  system  on  T3 ,  the  pathlength  between 
the  bending  magnets  will  be  different  from  the  central 
trajectory  and  the  electron  would  end  up  at  a  different 
location  downstream.  The  three  electron  trajectories 
illustrate  the  change  in  pathlength  due  to  the  effects  of 
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divergence  and  percent  momentum  change.  This  is  analogous 
to  the  change  in  the  longitudinal  extent  of  a  picobunch  as 
it  traverses  the  system. 

2 .   TRANSPORT  Equations 

The  TRANSPORT  equations  for  location  E  are  derived 
in  Appendix  B.   These  equations  are: 

xE(t)  =  (2p0tana+  Lq/cos2  a  +  L2)  9Q 

+  (-2P0  +  2  Pq/cos  a  +  L1sinc/cos2a)  6Q,     (M-13) 

lE(t)  =  (-2p0  +  2  pQ/cosa  +L1sina/cos2a)0o 

+  (-2p0a+  2pQtana  +  L1tan2a)50  +  10.     (M-15) 

It  can  be  shown  that  the  TRANSPORT  equations  reduce  to  the 
Buskirk  equations  letting  the  effective  pathlength  through 
the  bending  magnet,  L,  and  the  initial  one-half  longitudinal 
extent  of  the  picobunch,  10,  equal  zero.  The  variables  x 
and  1  correspond  to  dx  and  ds  respectively  in  Equations  2-1 
and  2-2.  TRANSPORT  analysis  requires  the  phase  space 
representation  of  the  input  beam.  For  a  normalized 
horizontal  beam  emittance  of  8  mm-mrad  (revised  MKIII  IRFEL 
parameter)  with  a  beam  energy  of  4  0  MeV,  assume  the  2-D 
phase  space  of  Figure  2-9.  The  points  a-d  corresponding  to 
the  extreme  values  of  the  phase  space  ellipse  are  to  be 
plotted.  This  will  yield  a  physical  representation  of  what 
the  detector  will  see  at  the  location  of  interest. 
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Figure  2-9   2-D  Phase  Space 

Substituting   the   system   parameters   with   10   =   .0005, 
Equations  M-13  and  M-15  become: 


X(t)  =  .5129680  +  .739766    +  xQ  , 


(2-3) 


l(t)  =  .O73980Q  +  .139406o   +  .0005  , 


(2-4) 


where  x0,  x(t)  ,  10,  l(t)  are  in  mm;  60  in  mrad,  6Q  in 
%  AP/P.  The  plot  of  x(t)  vs  l(t)  for  the  points  in  Figure 
2-9  is  shown  in  Figure  2-10.  Tabulated  values  are  shown  in 
Table  2-2.  The  center  ellipse  corresponds  to  a  mono- 
energetic  picobunch  with  divergence  in  the  xz  plane.  The 
divergence  alone  causes  a  significant  increase  in  the 
picobunch  extent.  A  picobunch  of  longitudinal  extent  1 
micron  grows  to  a  longitudinal  extent   of   approximately   14 
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microns  at  location  E.  The  right  side  ellipse  represents 
the  same  picobunch  having  one-tenth  of  one-percent  greater 
momentum  and  the  left  side  ellipse  represents  the  picobunch 
with  one-tenth  of  one-percent  less  momentum.  A  non- 
monoenergetic  picobunch  with  50  =  +/-  .1%  (revised  MKIII 
IRFEL  parameter)  would  have  a  total  longitudinal  extent  of 
=  46  microns  at  location  E.  A  coherent  TR  effect  would  not 
be  expected  under  these  conditions.  Even  in  the  mono- 
energetic  case,  the  beam  divergence  would  have  to  be 
greatly  reduced  to  yield  a  longitudinal  extent  of  <=  3 
microns.  The  TRANSPORT  system  equations  for  detector 
placement  midway  between  locations  B  and  C  are: 

xc(t)  =  (cosa)x0  +  (p0sina  +  L^/cosajSQ 

+  [p0(l-cosa)  +  L]_tana]6  Q  ,         (M-5) 

IcCt)  =  (-sim)x0  -  p0(l-cosa)60  +  1Q 

-  p0(a-sina)6  0.  (M-7) 

Substituting  the  known  parameters  with  L^  =  163.52  mm  these 
equations  become: 

x(t)  =  .98079xo  +.204826o  +  .36278  6q    ,  (2-5) 

l(t)  =  -.19509xo  -  .00375Go  -.002466o  +  .0005  .     (2-6) 
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A  plot  of  l(t)  vs  x(t)  for  the  phase  space  of  Figure  2-9  is 
shown  in  Figure  2-11.  Tabulated  values  are  shown  in  Table 
2-3. 

Inspection  of  Figure  2-11  reveals  a  shift  in  the 
orientation  of  the  ellipses  with  respect  to  Figure  2-10.  An 
11.25  degree  rotation  of  the  thin  foil  about  the  y-axis  will 
minimize  the  total  longitudinal  extent  of  the  original 
picobunch  as  seen  by  the  thin  foil.  The  effective 
longitudinal  extent  of  the  original  picobunch  is 
approximately  23  microns  for  constant  x(t).  Moving  the 
detector  further  upstream  to  a  location  10  mm  downstream  of 
the  first  bending  magnet  the  TRANSPORT  system  equations  are: 

X(t)  =  .98079xo  +  .048306o  +  .05742SQ  ,  (2-7) 

l(t)  =  -.19509xo  -  .003759o  -  .00246<5Q  +.0005  .     (2-8) 

A  plot  of  l(t)  vs  x(t)  for  this  detector  location  is  shown 
in  Figure  2-12.  Tabulated  values  are  shown  in  Table  2-4. 
Note  that  Equation  2-6  and  Equation  2-8  are  identical,  thus 
the  longitudinal  extent  of  the  picobunch  is  constant  between 
B  and  C.  However,  the  effective  picobunch  extent  will 
increase  between  B  and  C.  Figure  2-12  reveals  that  a 
momentum  spread  of  one-tenth  of  one-percent  will  not  cause 
significant  increase  in  the  effective  longitudinal  extent  of 
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the  picobunch.    The  non-monoenergetic  ellipses  appear  to 
closely  coincide  with  the  central  momentum  ellipse. 

It  is  apparent  from  equations  M-7  and  M-15  that  the 
change  in  the  longitudinal  extent  of  the  original  picobunch 
is  due  to  the  initial  horizontal  displacement  of  the 
electron,  x0,  the  divergence  of  the  electron,  60,  and  the 
momentum  increment,  5Q;  all  with  respect  to  the  assumed 
central  trajectory.  A  scaled  plot  is  used  to  illustrate  the 
relative  effect  of  each  factor  contributing  to  the  change  in 
the  horizontal  and  longitudinal  extent  of  the  picobunch  as 
it  traverses  the  system.  Substituting  the  known  parameters, 
the  TRANSPORT  equations  to  be  plotted  reduce  to  the 
equations  of  Table  2-1. 

TABLE  2-1 

EQUATIONS  FOR  TRAJECTORY  PARAMETERS 
AT  VARIOUS  LOCATIONS 

Detector  10  mm  from  bending  magnet  #1: 

X(t)  =  .98079xo  +  .O483O0O  +  .05742  60 

l(t)  =  -.19509xo  -  .OO3750O  +  1Q  -  .00246  6Q 

Detector  halfway  between  bending  magnets  #1  and  #2: 
x(t)  =  .98079xo  +  .204829o  +  .3627860 
l(t)  =  -.19509xo  -.0037500  +  10  -.00246So 

Detector  at  location  E: 

x(t)  =  xQ  +  .51296e0  +  .7397660 
l(t)  =  .07398eo  +  10  +  .139406o 


30 


Figure  2-13  is  an  unsealed  representation  of  the 
effect  of  the  initial  horizontal  displacement,  x0,  on  the 
longitudinal  extent  of  the  beam.  The  beam  parameters  are 
specified  by  the  given  2-D  phase  space.  Figure  2-13 
illustrates  that  the  initial  horizontal  displacement  with 
respect  to  the  central  trajectory  will  have  an  effect  on  the 
longitudinal  extent  of  the  beam  between  the  magnets  but  will 
have  no  effect  beyond  the  second  magnet  at  location  E. 

Figure  2-14  is  a  scaled  representation  of  the  effect 
of  the  divergence,  0O  on  x(t)  and  l(t).  The  effect  of  the 
divergence  term  on  l(t)  is  small  over  the  regions  between 
the  magnets  but  increases  beyond  the  second  bending  magnet. 
Note  that  the  change  in  the  longitudinal  extent  is  measured 
with  respect  to  the  beam  direction. 

Figure  2-15  illustrates  the  effect  of  the  momentum 
increment,  6Q/  on  x(t)  and  l(t). 

Comparing  Figures  2-14  and  2-15,  it  appears  that  the 
divergence  and  momentum  increment  have  roughly  the  same 
effect  in  debunching  the  beam  picobunches  as  they  traverse 
the  system.  The  previous  three  figures  and  Figures  2-10,  2- 
11,  and  2-12  lead  to  the  physical  representation  of  the 
beam  picobunches  shown  in  Figure  2-16. 


31 


w 


=—  < 


C 


c 
o 


4-1 

o 
u 

<D 
4-1 
4-1 

w 


I 

fN 
<D 

u 

•H 

fa 


32 


w 


4-> 

n 
c 
<a 

+j 

c 
o 

c 

CD 
4-1 

o 

+J 
u 

0) 

4-1 
m 
w 


U 


J3T 


33 


TDI 


w 


T3 

C 


c 
o 


o 
u 

w 


in 

I 

(N 

Q) 
M 

en 

•H 

Cm 


Hh- 


34 


4-> 

8 

u 

tn 

•H   -P     1 

> 

(U  S   q) 

en  en  q 

•H 

idiS 

U 

g  >  s 

o 

iff 

•H 
-M 

c  c  c 

0) 

sss 

en 
1 

KWII 

till  Uljl  III  llh 


x: 
u 
c 

13 
O 
U 
•H 
CU 

I 

ffl 
M-t 

o 

c 

o 

•H 

+J 

c 
a; 
w 

<D 
U 

Cu 
<D 


U 
•H 
W 
>i 

x: 


rH 

I 

H 

Cn 
•H 
Cm 


ffiffiHH] 


35 


TABLE  2-2 
BEAM  PARAEMTERS  FOR  DETECTOR  AT  LOCATION  E 

Point      fx.  B]      %Ap/p 

a         (1/0)        0 

.1 
-.1 

b         (0, .1)       0 

.1 
-.1 

c         (-1,0)       0 

.1 
-.1 

d         (0,-.l)      0 

.1 
-.1 


1ft)    mm 

x ( t )    mm 

.001 

1.000 

.015 

1.074 

-.013 

.926 

.008 

.051 

.022 

.125 

-.006 

-.023 

.001 

-1.000 

.015 

-.926 

-.013 

-1.074 

-.006 

-.051 

.008 

.023 

-.020 

-.125 

TABLE  2-3 

BEAM  PARAMETERS  FOR  DETECTOR  HALFWAY 
BETWEEN  BENDING  MAGNETS 


Point 


fx,  6). 

%AP/P 

1ft)    mm 

x(t)    mm 

(1,0) 

0 

-.194 

.981 

.1 

-.194 

1.017 

-.1 

-.194 

.945 

(0/.1) 

0 

.001 

.020 

.1 

0 

.057 

-.1 

.001 

-.016 

(-1,0) 

0 

.196 

-.981 

.1 

.196 

-.945 

-.1 

.196 

-1.017 

(0,-.l) 

0 

.001 

-.020 

.1 

.001 

.016 

-.  1 

.002 

-.057 
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TABLE  2-4 

BEAM  PARAMETERS  FOR  DETECTOR  10  mm 
DOWNSTREAM  OF  FIRST  BENDING  MAGNET 


Point 


(x,6) 

%AP/P 

1ft)    mm 

x ( t )    mm 

(1,0) 

0 

-.194 

.981 

.1 

-.194 

.987 

-.1 

-.194 

.975 

(0,-1) 

0 

.001 

.005 

.1 

0 

.011 

-.1 

.001 

-.001 

(-1,0) 

0 

.196 

-.981 

.1 

.196 

-.975 

-.1 

.196 

-.987 

(o,-.i) 

0 

.001 

-.005 

.1 

.001 

.001 
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III.   RESULTS  AND  CONCLUSIONS 

The  axially  modulated  beam  of  the  MKIII  IRFEL  was 
represented  as  disks  of  charge  (picobunches)  as  shown  in 
Figure  1-4.  If  these  picobunches  were  to  retain  their  axial 
modulation  and  pass  through  a  thin  foil,  coherent  optical 
transition  radiation  would  be  generated. 

The  analysis  of  the  TRANSPORT  system  equations  for  tr~ 
detector  locations  of  Figure  2-5  reveal  that  significant 
debunching  (of  more  than  3  microns  in  the  axial  direction) 
will  occur  at  these  detector  locations  due  to  the  effects  of 
a  non-monoenergetic,  divergent  beam. 

For  detector  placement  at  the  four-way  cross  location 
(location  E  of  Figure  2-5)  the  increase  in  the  picobunch 
extent  is  approximately  46  microns  (Figure  2-10) .  For  the 
detector  location  halfway  between  bending  magnets  #1  and  #2 
(Figure  2-5)  the  effective  picobunch  extent  is  reduced  to  2  3 
microns  (Figure  2-11)  with  the  detector  rotated  11.25 
degrees  about  the  y-axis.  Placing  the  detector  further 
upstream  has  the  overall  effect  of  reducing  the  effective 
picobunch  extent.  The  axial  modulation  of  the  electron  beam 
of  the  MKIII  IRFEL  can  be  measured  at  a  detector  location  10 
millimeters  downstream  of  the  first  bending  magnet  with  the 
detector  rotated  11.25  degrees  about  the  y-axis.  At  this 
location  there  will  be  a  2  micron  increase  in  the  effective 
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picobunch  extent  (Figure  2-12) .  In  theory,  detector 
placement  at  this  location  will  enable  the  desired  OTR 
effect  to  be  observed.  This  location,  however,  is  not 
easily  accessible  due  to  structural  hardware  upstream  of  the 
bellows  of  Figure  2-4.  In  addition,  the  bellows  location 
and  waveguide  material  would  have  to  be  modified  in  order  to 
mount  a  detector  of  the  type  illustrated  in  Figure  2-6.  To 
situate  a  detector  at  this  location  it  will  be  necessary  to 
make  hardware  modifications  to  accomodate  the  experimental 
arrangement. 

The  normalized  horizontal  and  vertical  emittance 
parameters  for  this  experimental  proposal  were  8  mm-mrad  and 
4mm-mrad  respectively.  The  momentum  spread  was  one-tenth  of 
one  percent.  A  reduction  in  these  parameters  will  yield  a 
reduced  picobunch  extent  at  any  detector  location. 

In  addition  to  the  Stanford  MKIII  IRFEL,  there  are 
several  IRFELS  that  are  potentially  suitable  for  an 
experiment  of  this  nature.  Potential  candidates  for  this 
experiment  include: 

(1)  The  Free  Electron  Laser  Driven  by  the  NBS  Microtron 

(2)  Stanford   Superconducting  Accelerator/Free   Electron 
Laser  SCA/FEL 

(3)  Lawrence   Livermore   National   Laboratory   ATA/TRW 
Wiggler. 

The   NBS   CW   Microtron   Driven   FEL   (currently   under 

development)  has  the  design  goals  of  a  normalized  transverse 

emittance  less  than  10  mm-mrad  and  a  momentum  spread  of  less 
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than  .03%  for  operation  at  0.2-3.0  microns.  This  momentum 
spread  is  one  order  of  magnitude  better  than  the  MKIII  IRFEL 
due  to  the  inherent  acceleration  process  of  a  microtron. 
The  hardware  and  bending  magnet  configuration  downstream  of 
the  wiggler  exit  need  to  be  investigated  and  modeled  in 
order  to  determine  a  suitable  detector  location  to  measure 
the  electron  beam  modulation. 

The  Stanford  Superconducting  Accelerator/Free  Electron 
Laser  (SCA/FEL)  has  a  normalized  transverse  emittance  of  5 
mm-mrad  (10  mm-mrad  with  a  high  brightness  injector)  and  a 
momentum  spread  of  approximately  .04%  for  a  fundamental 
wavelength  of  3u.  [Ref.8]  The  hardware  and  bending  magnet 
configuration  downstream  of  the  wiggler  exit  also  need  to  be 
investigated  and  modeled  in  order  to  determine  suitability 
for  a  proposed  experiment. 

The  Advanced  Test  Accelerator/TRW  Wiggler  FEL  utilizes 
induction  linac  technology  and  the  laser  amplifier  FEL 
configuration  to  amplify  a  carbon-dioxide  laser  at  10.6 
microns.  This  system  may  be  more  suitable  for  detector 
placement  due  to  the  absence  of  optical  cavity  mirrors  and  a 
portable  beam  dump. 

The  TRANSPORT  equations  were  generated  using  a  first 
order  approximation  and  an  assumed  beam  spot  size  of  radius 
1  mm  at  the  entrance  pole  face  of  the  first  bending  magnet. 
The  accuracy  of  the  equations  could  be  improved  by  utilizing 
the  second  order  fitting  capability  of  TRANSPORT. 
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The  system  of  Figure  B-3  was  modeled  starting  at 
location  A,  assuming  a  beam  spot  size  of  radius  1  mm.  A 
more  accurate  model  would  incorporate  the  effects  on  the 
electron  beam  envelope  within  the  wiggler.  TRANSPORT  does 
not  have  a  designated  "wiggler"  element  however,  the  wiggler 
can  be  modeled  using  drift  spaces  and  the  "Arbitrary 
Transformation"  (type  code  14)  input  option.  This  model 
should  start  at  the  wiggler  center  where  the  beam  waists  xw 
and  yw  are  known.   [Ref.  10] 
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APPENDIX  A 
BEAM  TRANSPORT  THEORY 

1.   PHASE  SPACE 

The  transport  of  a  charged  particle  beam  through  a 

system  of  magnets  and  drift  spaces  involves  the  study  of  a 

fundamental  conservation  law  known  as  Liouville's  theorem: 

.  .  .  under  the  action  of  forces  which  can  be  derived  from 
the  Hamiltonian,  the  motion  of  a  group  of  particles  is 
such  that  the  local  density  of  the  representative  points 
in  the  phase  space  remains  everywhere  constant.   [Ref.  11] 

Liouville's  theorem  advises  the  reader  that  he  must  look  at 
the  collective  group  of  particles  rather  than  an  individual 
particle  in  order  to  access  the  overall  nature  of  beam 
transport.  A  single  particle  moving  in  a  three  dimensional 
coordinate  system  is  completely  specified  if  its  position 
and  momentum  components  are  known,  namely  x,  y,  z,  Px,  Py, 
Pz .  It  is  convenient  to  express  this  information  in  a  six- 
dimensional  space  referred  to  as  phase  space.  A  beam 
comprises  all  individual  particles  within  this  phase  space, 
neglecting  electrostatic  repulsion  and  spin  dependent 
effects.  rhe  latter  two  effects  would  require  an 
unmanageable  phase  space  of  dimension  six  times  the  number 
of  particles  in  the  beam  plus  the  additional  dimensions 
corresponding  to  the  rotational  degrees  of  freedom  of  the 
individual  particles.  Any  system  of  conservative  forces 
acting  on  a  beam  possesses  a  Hamiltonian.   The  macroscopic 
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magnetic  forces  along  the  beam  trajectory  of  the  proposed 
experiment  are  conservative  thus  the  beam  must  obey 
Liouville's  theorem.  The  concept  of  a  six-dimensional  phase 
space  is  a  valid  and  useful  utility  in  the  study  of  beam 
transport  since  space  charge  effects  become  weaker  as  beam 
energy  increase.  The  six-dimensional  phase  space 
representing  all  particles  in  the  beam  is  considered  a 
hypervolume  that  can  change  shape  but  must  maintain  constant 
volume.  The  phase  space  analysis  of  the  beam  can  be  further 
simplified  by  the  choice  of  the  coordinate  axis.  If  the 
three  components  of  motion  are  mutually  independent  in  real 
space,  the  motion  is  confined  to  the  three  planes  (x,  Px)  , 
(Y/  py)  /  anc*  (z/  pz)  •  T^e  s^x  dimensional  phase  space 
hypervolume  has  in  effect  been  reduced  to  three  two- 
dimensional  phase  space  areas.  A  further  simplification  in 
the  analysis  is  possible  by  assuming  the  axial  momentum  of 
the  beam  is  constant.  This  is  a  valid  assumption  downstream 
of  the  accelerator  sections  of  the  beam  trajectory.  The 
angular  divergence  of  a  particle  relative  to  the  beam  axis 
is  equal  to  the  ratio  of  the  transverse  and  axial  momenta. 
Considering  the  transverse  x-direction  the  axial  divergence 
would  be  Px/Pz  or  dx/dz.  Dimensional  phase  space  (x,Px)  is 
now  represented  as  (x, dx/dz)  where  dx/dz  is  the  angular 
divergence,  a  directly  observable  quantity.  Similarly,  the 
transverse  y-direction  phase  space  representation  is 
(Y/dy/dz) .   The  longitudinal  phase  space  is  represented  as 
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(l/dPz/Pz)  where  1  denotes  the  longitudinal  beam  extent  and 
dPz/Pz  denotes  the  fractional  momentum  spread.  [Ref..  11] 

2 .   EMITTANCE 

The  quality  of  any  electron  beam  is  measured  in  terms 
of  its  emittance.  Consider  a  single  electron  propagating  in 
the  z-direction  as  shown  in  Figure  A-l. 


electron  having  off  axis 
component  of  velocity 


+  z 


electron  on  beam  axis 


Figure  A-l  Single  Electron  Motion  in  the  Horizontal  Plane 

Emittance  is  conveniently  expressed  in  terms  of  the  phase 
space  parameters  x  and  dx/dz.  Using  the  small  angle 
approximation  dx/dz  can  be  represented  as  0.  Electron  beams 
produced  by  accelerators  have  elliptical  phase  spaces  as 
shown  in  Figure  A-2  where: 

0max  =  the  maximum  divergence  of  an  individual  electron 
xmax  =  the  maximum  horizontal  extent  of  a  beam  electron. 
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Figure  A-2   Elliptical  Phase  Space 

In  accordance  with  Liouville's  theorem,  phase  space  is 
conserved.  The  electrons  are  within  the  perimeter  of  the 
ellipse.   Emittance  is  mathematically  defined: 


e  =  (area  of  phase  space ) /^    . 


(A-l) 


Normalized  emittance  is  defined: 


:n  =  ^e  > 


(A-2) 


where  y  is  the  ratio  of  the  particle  energy  to  the  rest  mass 
energy  and  3=1.  The  area  of  the  elliptical  phase  space  is 
irXmaxOmax  therefore: 
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^h  =  Xmax^max. 

For  example:    The  Stanford  MKIII  FEL  has  the  following 
normalized  parameters  [Ref.  12]: 

horizontal  phase  space  area  =  8  -nrnm-mrad  , 
vertical  phase  space  area  =  4  t  mm-mrad  . 

therefore  assuming  E  =  40MeV: 

£  nh  =  8  itum-mrad;  £  h  =  • 10  mm-mrad  , 
e  nv  =  4  mm-mrad;  e  v  =  .05  mm-mrad  . 

Assuming  elliptical  phase  spaces,  we  can  assign  approximate 
values  for  Xmax  and  Ymax  and  the  corresponding  values  for 
6max  and  <j>  max. 

Letting  Xmax  =  Ymax  =  1mm  and  energy  equal  to  4  0MeV  with 
B  =  1  the  corresponding  values  of  max  and  max  are 
determined  to  be: 

Xmax  =  1  mm;  0  max  =  .10  mrad  , 

Ymax  =  1  mm;  *  max  =  .05  mrad  . 
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3.   BEAM  WAISTS 

In  accordance  with  electron  beam  transport  design  [Ref. 
10],  the  electron  beam  is  normally  focussed  at  the  center  of 
the  wiggler  of  the  FEL.  This  yields  a  beam  waist  in  the 
center  of  the  wiggler  as  shown  in  Figure  A-3. 
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Figure  A-3   Beam  Focus  at  Center  of  Wiggler 


where 


xw  =  horizontal  beam  waist 

yw  =  vertical  beam  waist 

Ly,  =  length  of  wiggler  (1  meter) 

The  horizontal  beam  waist  size  is  related  to  the  length  of 
the  wiggler  and  the  horizontal  emittance  by  the  relation: 


xw  =   Lweh/2 


(A-3) 
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The  vertical  beam  waist  size  is  related  to  the  vertical 
emittance  and  the  intrinsic  betafunction,  3  w,  by  the 
relation: 


Yw  =  /£vew  (A-4) 

where 


w  =  electron  radius  of  curvature  at 
center  of  wiggler. 
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APPENDIX  B 

TRANSPORT-A  COMPUTER  PROGRAM  FOR  DESIGNING 
CHARGED  PARTICLE  BEAM  SYSTEMS 


1.   INTRODUCTION  TO  TRANSPORT 

The  study  of  static-magnetic  beam  transport  systems  led 
to  the  development  of  the  original  first-order  TRANSPORT 
computer  program.  TRANSPORT  was  written  in  BALGOL  at 
Stanford  Linear  Accelerator  Center  (SLAC)  in  1963.  It  has 
been  refined  in  the  succeeding  years  by  contributions  from 
various  laboratories  throughout  the  world.  Subsequent 
improvements  include  the  translation  into  FORTRAN  and  a 
second-order  fitting  capability.  The  version  of  TRANSPORT 
used  in  this  thesis  is  the  PC  version  dated  May,  1985  from 
D.  Carey,  Fermi  National  Accelerator  Laboratory  (FNAL) , 
compiled  with  the  IBM  version  of  professional  fortran  by 
Ryan-McFarland.   [Ref.  13] 

TRANSPORT  utilizes  matrix  theory  methods  to  step 
sequentially  through  a  beam  line  system  consisting  of 
magnets,  drift  spaces,  and  special  magnet  configurations. 
The  beam  line  system  and  the  input  beam  parameters  are 
specified  by  the  user  and  an  output  beam  matrix  and  overall 
system  transfer  matrix  are  generated.  Each  individual  step 
through  the  system  is  concidered  an  element.  Each  element 
is  given  by  a  sequence  of  items  separated  by  spaces  and 
terminated  by  a  semicolon.    These  items  are  a  type  code 
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number,  a  vary  field,  the  physical  parameter  and  optional 
labels.  The  type  code  number  identifies  the  element  (i.e., 
magnet,  drift  space,  pole  face  rotation)  .  A  vary  field 
allows  the  user  to  adjust  the  physical  parameters  of  an 
element  if  the  problem  involves  fitting.  Physical 
parameters  are  quantities  which  describe  the  physical 
element.  Optional  labels  assist  in  identifying  beam 
elements  and  parameters  in  the  output.  A  summary  of 
TRANSPORT  type  codes  is  shown  on  page  77.   [Ref.  13] 

The  initial  objective  in  writing  the  input  program  will 
be  to  specify  the  beam  line,  establishing  the  central 
trajectory.  The  system  to  model  by  TRANSPORT  is  shown  in 
Figure  B-l.  This  system  consists  of  two  bending  magnets  and 
five  drift  lengths  where 

L0  =  initial  drift  length  , 

L  =  effective  pathlength  through  bending  magnet  , 

Lx  =  drift  length  1  , 

L2  =  drift  length  2  , 

a   =  bending  angle  . 

2.   SAMPLE  PROGRAM  DEVELOPMENT  [Ref.  13] 

Figure  B-2  is  a  sample  input  program  that  establishes 
the  beam  line  of  Figure  B-l,  assuming  the  initial  drift 
space,  L0  equals  zero.  Refer  to  Figure  B-2  for  the 
remainder  of  this  section. 
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Figure  B-l   System  Configuration 
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mm ' 

0. 

1 

MeV 

0 

001 

mm ' 

0. 

001 

mm' 

0 

.  1 

(1)  'STANFORD  MKI II  IRFEL:  REVISED  10/16/87' 
0 

(2)  (CARTESIAN  COORDINATE  SYSTEM;  x.y.z) 

13.      'COOR'      12.   ; 

(3)  (SET  UNITS/scaled  units) 
15.      ' xymm'      1 . 
15.      '  EM*»V '     11. 
15.       'lom;'      8. 
lb.      ' leng'      5 . 

(4)  (BEAM  PARAMETERS  AT  FIRST  BENDING  MAGNET  POLE-FACE) 
(NORMALIZED  BEAM  EMITTANCE:  Enh=8min-mrad .  Env=4mm-mrad ) 
1.0    'BEAM'   1.0    0.1    1.0    0.05    .0005   0.0    40.0 
( monoenergetic  beam) 

(5)  1G.      5.    5.0    ;   (half -gap), 
(default  value  of  Kl=.5) 

(6)  (FIRST  BENDING  MAGNET) 
?..  0 .   ; 

4.     'BNP1'     38.347      6.33175     0.   ; 

2.  11 .25    ; 

13.    1  .   ;  (PRINT  rr.;AM  MATRIX ) 
13.    4.   ;  (PRINT  XFORM  MATRIX) 

(7)  (DRIFT  SPACE  til  ) 

3.  'DFT1'    327.03   ; 
13.    1 .   ; 

13.    4.   ; 

(8)  (SECOND  BENDING  MAGNET) 

20.     180.     'LEFT'  ;   (st.d  coord  rot  w/11.25  bend) 

6.83175    0.   ; 


(9) 


2.            11.25 

; 

4.            ' BNDZ' 

38. 347 

2 .           0 .       ; 

20.            -130. 

\ 

13.            1.       ; 

13.            4 .       ; 

(DRIFT    SPACE 

«2) 

3.            'DFT2* 

95.30 

13.            1 .       ; 

13.           4 .       ; 

(10)     SENTINEL 


Figure  B-2   Sample  Input  Program 
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The  individual  elements  and  there  associated  type  codes  will 
be  explained  in  the  succeeding  paragraphs. 

(1)  'STANFORD  MKIII  FEL:  revised  10/16/87'  is  the 
arbitrary  title  of  the  sample  input  program.  The  zero 
immediately  below  the  title  indicates  the  start  of  a  new 
problem. 

(2)  A  cartesian  coordinate  system  was  chosen  for  the 
proposed  experiment  design  problem.  The  coordinate  system 
is  established  using  the  type  code  13.  The  statement  [13. 
12.  ;]  will  yield  a  beam  line  reference  trajectory  that  will 
start  at  the  origin  and  proceed  along  the  positive  z-axis. 
The  y-axis  will  point  up  and  the  x-axis  will  point  to  the 
left.  Statements  enclosed  by  (  )  are  comment  statements  to 
assist  in  following  the  code. 

(3)  The  scaling  of  units  is  accomplished  via  type  code 
15.  All  lengths  in  Figure  B-l  are  to  be  specified  in 
millimeters  and  energy  is  to  be  specified  in  MeV.  TRANSPORT 
has  the  standard  set  of  units  as  shown  on  page  78. 
Conversion  factors  are  applied  to  code  digits  1,  11,  8,  and 
5  to  establish  an  appropriate  unit  system.  A  summary  of 
TRANSPORT  type  codes  is  found  on  page  77  . 

(4)  The  input  beam,  denoted  by  the  type  code  1,  is 
specified  in  terms  of  its  phase  space  and  average  momentum 
parameters.  The  phase  space  input  is  given  in  terms  of  the 
semi-axis  of  a  six-dimensional  beam  ellipsoid  representing 
the  phase   space  variables  x,   0,   y,  <j>  ,   1,  and  6  where 
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x  =   one-half  the  horizontal  beam  extent  (1.0mm)  , 
8  =   one-half  the  horizontal  beam  divergence  (O.lmrad)  , 
y  =   one-half  the  vertical  beam  extent  (1.0mm)  , 
<J>  =   one-half  the  verticle  beam  divergence  (.05mrad)  , 
1  =   one  half  the  longitudinal  beam  extent  (.0005mm)  , 
6  =   one-half  the  momentum  spread  (0.0%)  , 
P  =   average  momentum  of  the  beam  (4  0MeV/c) . 
NOTE:   The  above  units  were  scaled  using  type  code  15. 
(5)  and  (6)   Bending  magnet  parameters  are  specified  in 
terms  of  the  wedge  magnet,  type  code  4;  pole-face  rotation, 
type  code  2;  and  a  gap  setting,  type  code  16.    A  wedge 
bending  magnet  implies  that  the  central  trajectory  of  the 
beam  enters  and  exits  perpendicularly  to  the  pole  face 
boundaries.   A  wedge  magnet  can  be  modified  into  a  square 
edge  magnet  using  the  pole  face  rotation  element  on  either 
the  entrance   or  exit   face  of  the  magnet.    The  field 
boundaries  for  bending  magnets  are  shown  on  page  73.    To 
modify  a  wedge  bending  magnet  into  a  square  for  a  bend  to 
the  right,    B]_  is  set  to  zero  and  32  is  set  to  11.25 
degrees.     The   pole   face   rotation   type   code   2.   must 
immediately  precede  or  immediately  follow  the  bending  magnet 
(type  code  4.0).   There  are  four  first  order  parameters  to 
be  specified  for  the  wedge  magnet;   the  type  code,   the 
effective  length  of  the  central  trajectory,   the  central 
field  strength  and  the  field  gradient.   The  bending  magnet 
equations  and  parameters  for  the  system  are  listed  below: 


54 


L  =  p0a  where 

L  =  effective  pathlength  through  magnet  , 

p0  =  bending  radius;  (19  5.30  mm)  , 

a  =  bend  angle;  (11.25  degrees)  , 

L  =  (195.30mm) (11. 25deg) (  rad/180deg)  =  38.347  mm  . 

(NOTE:  To  simulate  a  gradual  variation  of  the  fields  at 
the  ends  of  the  bending  magnets  a  fringing  field 
element  is  provided.) 

B(0)  =  33.356(P/p0)  where 

B(0)  =  central  field  strength  , 

P  =  central  momentum;  (40  MeV/c)  , 

B(0)  =  33.356(40MeV/c)/195.30mm  =  6.83175  kG  , 

n   =  magnetic  field  gradient  (n  =  0  assumed) . 

Code  [13.  1.  ;]  yields  the  printing  of  the  current  beam 

(a)  matrix.    Code  [13.  4.   ;]  yields  the  printing  of  the 

current  transformation  matrix,  Rl.    These  codes  will  be 

addressed  in  the  matrix  formulation  of  TRANSPORT. 

(7)  and  (9)   A  drift  length  (type  code  3.0)  is  a  field 

free  region  through  which  the  beam  passes.   It  is  denoted  by 

the  type  code  and  the  effective  drift  length: 

L]^  =  drift  length  1   (327.03  mm)  , 

L2  =   drift  length  2   (95.30  mm) 

(8)   The  second  bending  magnet  differs  from  the  first 

in  that  B-^  is  set  to  11.25  degrees  and  ffe    is  set  to  zero 

degrees.   In  addition,  a  coordinate  rotation  element  (type 

code  2  0.)  is  required  to  rotate  the  transverse  coordinates  x 

and  y  through  an  angle  about  the  z-axis.   A  bend  to  the  left 
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(looking  in  the  direction  of  beam  travel)  is  accomplished  by 
rotating  the  x,y  coordinates  by  18  0  degrees. 

(10)   The  word  "sentinel"  denotes  the  end  of  the  input 
program. 


3.   MATRIX  FORMULATION  OF  TRANSPORT  [Ref.  13] 
A.   Phase  Space  Vector 

At  any  specified  position  in  the  system  an  arbitrary 
charged  particle  is  represented  by  a  vector  X,  whose 
components  are  the  positions,  angles,  and  momentum  of  the 
particle  with  repect  to  the  central  trajectory.  The  X 
vector  is  represented: 

x 


X  = 


y 

i 

6 


where 
x 


=  the  horizontal  displacement  of  the  arbitrary  ray 
with  respect  to  the  assumed  central  trajectory  , 

G  =  the  angle  this  ray  makes  in  the  horizontal  plane 
with  respect  to  the  assumed  central  trajectory  , 

y  =  the  vertical  displacement  of  the  ray  with  respect  to 
the  assumed  central  trajectory  , 

$ =  the  vertical  angle  of  the  ray  with  respect  to  the 
central  trajectory  , 

1  =  the  path  length  difference  between  the  arbitrary  ray 
and  the  central  trajectory  , 

6  =  Ap/P;  the  fractional  momentum  deviation  of  the  ray 
from  the  assumed  central  trajectory. 
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B.   The  Transfer  Matrix  R 

The  passage  of  a  charged  particle  through  an  element 
(magnet,  drift  space,  rotation,  etc.)  may  be  represented  by 
Equation  B-l: 

X(l)  =  RX(O)  ,  (B-l) 

where 

X(0)  is  the  initial  coordinate  vector, 

X(l)  is  the  final  coordinate  vector, 

R    is  the  transformation  matrix  of  the  element. 

In  a  system  of  several  magnets,  pole  face  rotations,  and 
drift  spaces,  R  is  replaced  by  the  product  matrix  R(t) : 

R(t)  =  R(n) . . .R(3)R(2)R(1)  .  (B-2) 

R(t)  is  the  product  of  the  individual  element  matrices  of 
the  system.  Equation  B-3  represents  the  matrix  equation  for 
a  system: 

6 
Xi(t)  =  I     RijXj(O)  ,  (B-3) 

j=l 

where 

Xi  =  x, 

x2  =  e, 
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*3 

x4 


=  y. 


Xc  = 


x*  = 


1 ,  and 

6. 


This  product  matrix  is  automatically  calculated  by  the 
program  and  is  called  TRANSF0RM1.  TRANSF0RM1  can  be  printed 
when  desired  using  the  statement  [13.  4.  ;]  (e.g.,  paragraph 
(6)  )  .  For  static  magnetic  systems  possessing  midplane 
symmetry,  the  six  simultaneous  linear  equations  represented 
by  Equation  B-3  may  be  expanded  in  matrix  form  as  [Ref.  14]: 


x(t) 

e(t) 

y(t) 

<J>(t) 

i(t) 

6(t) 

_ 

Rll  R12  0  0  0 

R21  R22  °  °  ° 

0  0  R33  R34  0 

0  0 

R51  R52  0  0  1 

0  0  0  0  0 


R43   R44   ° 


R16 

R2  6 
0 

0 

R56 

1 


xo 

So 

Yo 

$o 

lo 

LV 

The  transformation  is  from  an  initial  position  t  =  0 
to  a  final  position  t  =  t  measured  along  the  assumed  central 
trajectory.  The  final  coordinate  vector  X(t)  is  the  matrix 
product  of  the  individual  coordinate  vector  and  the  overall 
system  transformation  matrix.  The  individual  R-^j  elements 
of  the  matrix  are  dimensionally  consistent  with  the  input 
units  specified  in  the  transport  code.  The  individual 
transformation  matrices  needed  to  represent  the  system  of 
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Figure  B-l  are  shown  on  pages  74-76.  It  is  desired  to 
determine  and  plot  x(t)  vs  l(t)  for  the  overall  system. 
Since  the  system  is  symmetric,  x(t)  and  l(t)  have  no 
dependence  upon  y0  or  fyQ,  thus  the  general  form  of  the 
tranformation  matrix  is  reduced  to: 


x(t) 
6(t) 
l(t) 
5(t) 


Rn  R12  ° 


R21   R22   0 


R51   R52   1 


R16 

xo 

R2  6 

8o 

R56 

lo 

1 

V 

The  wedge  bending  magnet  is  reduced  to: 


cos  a 

PoSin  a 

0 

P0(l-cosa) 

-hsin  a 

cos  a 

0 

sin  a 

-sin  a 

-p0(l-cos  a) 

1 

-  p0(  ct-sina) 

0 

0 

0 

1 

The  pole  face  rotation  matrix  is  reduced  to 


1 

0 

0 

0 

(tana)/Po 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 
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The  drift  space  matrix  is  reduced  to 


1 
o 

0 
0 


L 
1 
0 
0 


0 
0 

1 

0 


0 
0 
0 

1 


It  is  desired  to  determine  the  X  matrix  at  various  locations 
along  the  central  trajectory.  Figure  B-3  denotes  the  system 
locations  to  be  evaluated: 


Figure  B-3   System  Configuration  Evaluation  Points 


60 


The  matrix  representation  of  XB(t)  is  the  product  of  the 
pole  face  rotation  matrix,  wedge  bending  magnet  matrix,  and 
the  initial  coordinate  vector: 


xB(t) 

eB(t) 
IbW 


1  0   0   0 

tana/p0  10   0 

0  0   10 

0  0   0   1 


cos  a   pQsina   0   p0(l-cosa) 
-hsina   cos  a   0      sinct 
-sin  a~p0(l-cosa)  1  -p0(a-sina)j 
0  0      0'       1 


(pole  face  rotation) 


The  matrix  product  reduces  to: 


x 


Oj 


«o- 


(wedge  bending  magnet) 


xB(t) 

8B(t) 
lB(t) 
6B(t) 


cos  a 

p0sin  a 

0 

p0  (1-cos  a) 

0 

1/cos  a 

0 

tan  a 

sin  a 

"  Pc 

,(l-cosa) 

1 

-  Pq(  a-sin  a) 

0 

0 

0 

1 

lo 
9o 


The  matrix  representation  of  Xc(t)  is  XB(t)  multiplied  by 
the  first  drift  space  matrix: 


xc(t) 
6c(t) 

ic(t) 

6c(t) 


1      Lx  0       0 

0       1  0       0 

0       0  10 

0       0  0       1 

(first  drift) 


cos  a      f^sina        0  f^(l-cosa) 
0              1/cos  a        0  tan  a 

•sim  -pQ(l-cosa)    1  -p0(a-sina) 
0                     0                  0  1 

XB(t) 


x 
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The  matrix  product  reduces  to: 


xc(t) 
6c(t) 
lC(t) 

5c(t) 


cos  a  p0sina+L]_/cosa 
0  1/cosa 

■sina  -p0(l-cosa) 
0  0 


0  pQCl-cosaJ+LQ^tan   a 

0  tan  a 

1  -p0(a-sina) 
0         1 


x 


The  matrix  representation  of  XD(t)  includes  a  coordinate 
rotation  which  will  in  effect  bend  the  beam  to  the  left.  In 
matrix  form  this  is  represented  by  changing  the  bending 
radius  and  bending  angle  to  its  negative  value  (i.e.,  p0  -» 
-Pc  a  — =*~a)  •  The  wedge  bending  matrix  for  a  bend  to  the  left 
is  modified  as  shown  : 


cos  a      P0sin  a 


-hsina     cos  a 
sina     Po(l-cosa) 
0  0 


0 

p0  (1-cosa) 

0 

-sina 

1 

-p0(a-sina) 

0 

1 

The  pole  face  rotation  matrix  retains  its  identical  form 
because  tan(-a)  =  -tan ( a) •  Multiplying  Xc(t)  by  the  pole 
face  rotation  matrix  and  the  wedge  bending  magnet  matrix 
yields: 
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xD(t) 

1 

2  pQtancd-Lx/cos2  a 

0 

Zl 

9D(t) 

0 

1 

0 

0 

loCt) 

0 

Z2 

1 

Z3 

5D(t) 

0 

0 

0 

1 
J 

where: 

X, 


o 


Zl   =  j-2  p0   +    2pQ/cosa  +    L1sina/cos2a    , 

Z2    =    Zl    , 

Z3  =  -2p0a  +  2p0tana  +  I^tan2  a  . 

The  matrix  representation  of  XE(t)  is  XD(t)  multiplied  by 
the  second  drift  space  matrix: 


xE(t) 

9E(t) 
lE(t) 
6E(t) 


1 
0 
0 
0 


Z4 
1 

Z6 
0 


0 
0 

1 

0 


Z5 

0 

Z7 

1 


x 


where 


Z4    =    2  P0tan  a  +    L1/cos2  a  +    L2    , 


Z5   =   -2  p0   +    2  P0/cos  a  +    I^sin  o/cos2  a  , 


Z6    =    Z5    , 


Z7    =   -2  pQ  a  +    2  P0tan  a  +    I^tan' 


a 
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In  equation  form,  the  matrix  products  are  represented  as: 


xB(t) 

= 

(cosa)x0   +    (pQsina)  Qq   +  p    (1-cosa)  60 

(M-l) 

ee(t) 

= 

(1/cosa)  G0   +    (tana)  6q 

(M-2) 

iBCt) 

= 

(-sina)x0   -     ^(l-cosajBo  +   10   -    Pq  (  a  -   sina)Sc 

(M-3) 

6B  (t) 

= 

60 

(M-4) 

xc(t) 

= 

(cosa)x0  +    (Posina+   L^cos  a)  Qq  +    [p0(  1-cosa) 
+   Lxtan  a]  6q 

(M-5) 

€fc(t) 

= 

(1/cosa)  60  +    (tana)  5Q 

(M-6) 

lc(t) 

= 

(-sina)x0   -    p0(  1-cosa)  0O  +   lo   -    P0(a-sina)  6Q 

(M-7) 

6c(t) 

= 

«o 

(M-8) 

xc(t) 

= 

XB(t)     +    I^e^t) 

(M-5A) 

ec(t) 

= 

%(t) 

(M-6A) 

lC(t) 

= 

lB(t) 

(M-7A) 

5C(t) 

= 

6B(t) 

(M-8A) 

xD(t) 

= 

xo   +    (2  ftDtan  a+    Lq^/cos2  a)  %   +    (-2  Pq   +    2  pQ/coSx 
+    I^sin  c/cos2  a)  6q 

(M-9) 

eD(t) 

= 

6o 

(M-10) 

lD(t) 

= 

(-2  fo   +    2  to/cos  a+    I^sin  a/cos2  a)6  0   +    10 

+    (-2  poa  +    2p0tana  +    L1tan2a)6Q 

(M-ll) 

«D(t) 

= 

60 

(M-12) 
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xD(t 

0D(t 

6n(t 


=     (l/cosa)xc(t)     +     (posin    )  ^(t)     ~     Po(l-cosa)5c(t) 

(M-9A) 


=    (cosa)ec(t)    -    (siria)5c(t) 

=    (tana)xc(t)    +     pq  (l-cosa)e  c(t)    +    lc 
-   Po(a  -   sin    )  6c(t) 

=    6c(t) 


(M-10A) 
(M-11A) 

(M-12A) 


xE(t 
lE(t 

xE(t 

eE(t 

lE(t 
6E(t 


x0   +    (2p0tana  +   L^cos^    +   L2)e  0 

+    (-2p0   +    2p0/cosa+    Lxsina/cos2^  )  <5Q 


=    9o 

=  (-2p0  +  2  (^/cos  a  +  I^siro/cos2  a)  8   +  1Q 
+  (-2p0a  +  2p0tana  +L1tan2a)oQ 

=  6o 

=  xD(t)  +  L2eD(t) 
=  %(t) 

=  lQ(t) 

-  5D(t) 


(M-13) 

(M-14) 
(M-15) 

(M-16) 

(M-13A) 
(M-14A) 
(M-15A) 
(M-16A) 


C.   Beam  Matrix 

The  phase  ellipse  beam  (sigma)  matrix  is  used  to 
inform  the  user  of  the  beam  parameters  at  any  location  in 
the  system.  The  beam  matrix  will  represent  the  bundle  of 
particles  or  rays  constituting  the  phase  space  of  the  input 
beam  at  desired  locations  in  the  system.  Particles  of  the 
beam  are  assumed  to  lie  within  the  boundaries  of  the 
ellipsoid  with  each  point  within  the  ellipsoid  representing 
a  possible  ray.  The  equation  of  an  n-dimensional  ellipsoid 
may  be  written  in  matrix  form  as: 
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X(0)Ta (0)~1X(0)  =  1  ,  (B-4) 

where: 

X(0)T  =  transpose  of  coordinate  vector  X(0), 

a(0)   =  a  real  positive  definite,  symmetric  matrix. 

As  a  particle  passes  through  the  system,  it 
undergoes  the  matrix  transformation  of  Equation  B-l. 
Combining  this  transformation  with  the  equation  of  the 
initial  ellipsoid  and  using  RR_1  =  I  : 

X(0)T(RTRT"1)  a(0)"1(R"1R)X(0)  =  1  .       (B-5) 

Equation  B-5  can  be  rewritten  as: 

[RX(O) ]T[Ra(0)RT]"1[RX(0) ]  =  1  .  (B-6) 

The  equation  of  the  new  ellipsoid  after  the  transformation 
becomes: 

X(l)Ta(i)-lx(l)  =  1  ,  (B-7) 


where: 


a  (1)  =  RG(0)RT   .  (B-8) 
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Thus  the  beam  may  be  represented  at  any  point  in  the  system 
by  Equation  B-8.  Note:  a  (0)  is  the  input  beam  phase  space 
specified  by  the  user.  The  beam  matrix  is  printed  using  the 
statement  [13.  1.  ;].  The  projection  of  the  semi-axis  of 
the  ellipsoid  upon  each  of  its  six  coordinate  axes  is 
printed  in  a  vertical  array,  and  the  correlations  among 
these  components  are  printed  in  a  triangular  array.  The 
phase  ellipse  beam  matrix  in  the  printed  output  form  is 
shown  in  Figure  B-4. 
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where; 


r(ij)  --?iiii 


(p(ii)c(jj)] 


h 


As  a  result  of  the  fact  that  the  0  matrix  is  positive  definite,  the 
r(ij)  satisfy  the  relation 

|r(ij)|  £  1  . 


Figure  B-4  Beam  matrix 
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(Note:  units  will  be  labeled  as  specified  by  type  code  15.) 
The  individual  sigma  elements  are  defined: 


v'o     (11)  =  X 


max 


=  the  maximum  half-width  of  the  beam 
envelope  in  the  x (bend) -plane  at  the 
point  of  the  printout, 


So    (22)  =  Smax  =  tne  maximum  half -angular  divergence  of 

the  beam  envelope  in  the  x(bend)  plane, 

So    (33)  =  ymax  =  the  maximum  half -height   of   the   beam 

envelope, 

,/5"  (44)  =  cf^ax  =  the  maximum  half-angular  divergence  of 

the   beam   envelope   in   the   y (non-bend) 
plane, 

So    (55)  =  lmax  =  one-half  the  longitudinal  extent  of  the 

bunch  of  particles, 


So    (66)  =  6 


=  the  half-width  1/2 ( Ap/p)  of  the  momentum 
interval  being  transmitted  by  the  system. 


Note:  the  units  appearing  next  to  the  So  (ii)  in  the 
TRANSPORT  printout  are  the  units  chosen  for  coordinates  x, 
§  Yi    <$>  i    If  an<3  5  respectively. 


An  example  of  the  (x,6)  plane  ellipse  is  illustrated 
in  Figure  B-5. 

Consider  a  two  dimensional  (x,  0  )  plane  projection 
of  the  general  6-dimensional  ellipsoid.  Let: 


a   - 


^11 


°21 


a21 


a22 
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o       =  Ja 
max      v    22 


?int  =  ^22^-r21» 


CENTROID 


Q        r21  Xmax 


»-^ 


I  slope  = 


1       l'J22 


slope  =  r       /— — 
21    max 


*-  X 


11        max 


i 1         21         mt 


Figure  B-5  A  Two  Dimensional  Beam  Phase  Ellipse 

be  a  real,  positive  definite,  symmetric  matrix.   The  inverse 
of  the  a -matrix  exists  and  is  defined  by: 


a"1  =  1/c2 


°22 


-<2l 


-o 


21 


11 


where  e2  is  the  determinant  of 

The   two-dimensional   coordinate   vector   and 
transpose  are  represented  by: 

x 


its 


X  = 


e 


;   XT  =  (x,  e)  . 
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The  expansion  of  the  matrix  equation 


XTa"1X  =  1 


is  the  equation  of  an  ellipse,  therefore 


a  22x2  "  2a2ix9  +  ai:L  G2  =  e2  =  det 


The  area  of  the  ellipse  is  given  by: 

1/2 
A  =tt  (det  a  )  =  7rxmaxeint  =  TTXintemax  =  ve. 


The  correlation  between  x  and  8  (the  orientation  of  the 
ellipse)  depends  on  the  off  diagonal  term  a21-  This 
correlation  is  defined  as: 

a21 
r21  -  r12  " 


/alla22 

The  correlation,  r,  measures  the  tilt  of  the  ellipse  and  the 

intersection  of  the  ellipse  with  the  coordinate  axis.  The 

beam  ellipse  in  the  (x,  6  )  plane  will  be  valuable  in  the 

analysis  of  the  beam  divergence  at  various  points  in  the 
system  and  also  in  the  fact  that  it  reveals  the  change  in 
the  longitudinal  extent  of  the  beam  via  element  /o   (55) . 
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4.   INSTALLATION  OF  TRANSPORT  ON  THE  IBM  PC  [Ref.  15] 

The  following  requirements  and  procedures  are  necessary 
for  the  loading  and  running  of  TRANSPORT: 
Requirements : 

-  config.sys  file  must  have  files  >=  15;  buffers  >=  15. 

-  disk  set  A  containing  .EXE  and  .BAT  files  in  DOS/BACKUP 

-  64 OK  memory 

-  hard  disk;  lOM-byte  or  larger 
Procedures: 

-  move  to  the  C  directory:   C: [ENTER] 

-  create  a  TEMPO  subdirectory:  MD\TEMPO[ ENTER] 

-  move  to  the  TEMPO  subdirectory:   CD\TEMPO[ ENTER] 

-  RESTORE  .EXE  and  .BAT  files:  RESTORE  A:  C: [ENTER] 
(insert  disk  #1  and  #2  as  directed) 

-  move  to  the  C  directory:   CD\C: [ENTER] 

-  create  a  TRANS  subdirectory:   MD\TRANS [ENTER] 

-  move  to  the  TEMPO  subdirectory:   CD\TEMPO[ ENTER] 

-  copy  files  from  TEMPO  to  TRANS:   copy  *.*  \TRANS [ENTER] 

-  delete  files  from  TEMPO  subdirectory:   del  *.*[ ENTER] 

-  move  to  the  C  directory:   CD\C: [ENTER] 

-  remove  TEMPO  subdirectory:   RD\TEMPO[ ENTER] 
Limitations: 

-  the  PC  version  of  TRANSPORT  will  allow  for  1000  elements 
and  4500  data  values 

-  accuracy  is  four  significant  figures 
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To  run  TRANSPORT: 

-  TRANSGO.BAT  is  the  command  file  to  run  TRANS.EXE 

-  type  the  following: 

TRANSGO  FILE1  FILE2  FILE3 
FILE1:  input  data  file  name  (users  file  name) 
FILE2 :  output  print  file  name  (users  file  name) 
FILE3 :  output  summary  file  name  (users  file  name) 

-  example:   TRANSGO  MK3FELA  MK3FELA.0UT  MK3FELA.JNK[ ENTER] 

*  Note:  TRANSPORT  is  currently  loaded  in  the  IBM  PCs  in  the 
Physics  Library  and  Physics  Department  Office  under 
the  subdirectory  C:\TRANS.  X.  Maruyama  currently 
holds  the  TRANSPORT  software  on  5-1/4  disks. 
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Table  I:    Summary  of  TRANSPORT  type  codes 


PHYSICAL  BIB-CUT 

TYPE 

UJIE 

2nd 
ENTRY 

3rd 
ENTRY 

4th 
ENTRY 

5th 

ENTRY 

6th 
t>TRY 

7th 
EVi  RY 

8th 
LNTRY 

9  th 
ENTRY 

BEAM 

l.wwwO 

x(ca) 

e(mr) 

y(o») 

t(mr) 

I  (en) 

i  (percent) 

P. 

r.m.s.   AHOITION  TO 
BEAM  L\N\ELOFE 

1 .  wwwOO 

ox(os) 

69  (mr) 

6y(aa) 

6»(mr) 

6t(aaJ 

Ao(percent) 

6P 
(CeV/c) 

0 

IX)1£  FACE   ROTATION 

2.v 

rW.IF.DF 

H/IAIItlN  (degrees) 

1*1  FT 

3.v 

LENGTH    [metres) 

1EMDINC  MAGNET 

4.WV 

LENGTH   (metres) 

FIELD  (kG) 

FILLD  GRADIENT 
(n-value) 

qUADRUPOLE 

S.wO 

LENCTH  (netres) 

FIELD  (kG) 

HALF-APERTURE   (en) 

TRANSFORM  1   KTlATE 

6.0 

0.0 

1.0 

ntAKSHVM  2  irriMn 

6.0 

0.0 

2.0 

BEAM  (UflmiD  shut 

7. WWW 

SHIFT  (x)(cm) 

•3IIFT    (8)    (mr) 

SHIFT  (y)(cm) 

SHIFT    (»l(.r) 

SHIFT   (t)(cn) 

SHUT    («pcrcrnt) 

ALIOfrtNT  TOI-ERANCr 

8.  wwwO 

UISPLACLMTNT   (x)(cml 

ROTATION   (9l(mr) 

DISPUUJ1NT  (y)(cm 

ROTATION    (e)(mr) 

DISPLACEMENT    (l)(anl    ROTATION   (o)(mr) 

UIIC 

N1MUR 

REPEAT  CONTROL 

9.0 

N1KBER  OF 
REI'EAIS 

FITTING  QTBTRAIWrS 

10.0 

ll 

J 

HFSIPy  n  VAIUE  OF 
(I.J)   MATRIX 
IIJX  NIS 

ACCURACY  OF 
FIT 

Note:       •!    ll  used  f 

or    fitting 

■  beam  (o)  matrix  element.           *I    is  used   for  fitting  an  Rl  matrix  element. 

-   (1   ♦  20)   i 

s  used  for 

fitting  an  R2  matrix  element. 

ACCELERATOR 

11.0 

LENGTH  (metres) 

E  (energy  gain) 

(CeV) 

a   (phase   lag) 
(degrees) 

fWAVELENCTIO  (oa) 

BEAM 

(Rotated  Ellipse) 

12.0 

TIE   FIFTEEN  CDRRELAT 

IONS  AMONG  HOI  SIX  El 

EMEHTS  fi"hl»  entry  ssi 

st  be  preceded  by 

a  type  code  1.0  entry 

•  ) 

isTwotfrnr 
options 

13.0 

OlVTROL 
COUE  NUMBER 

ARBITRARY   R  MATRIX 

U.wvmO 

R(J.l) 

R(J,2) 

R(J,3) 

R(J,«) 

"(J. 5) 

R(J,6) 

J 

WITS  CONTROL 
( Transport 
lumens  ions) 

is.o 

CODE 

UNIT  SYMBOL 

SCALE  FACTOR 
(if   required) 

QUADRATIC  TERM 
Oh   RENUING  T1CLD 

16.0v 

1.0 

t(D   -   «<^)' 

c.  in  units  of  rransv 

orse   length   (cw) 

MASS  OF  PARTICLES 
IN  BEAM 

16.0 

3.0 

M/m  (dimensionless) 

m  -  mass  of  electron 

HALF-ArrRTURE  OF 
BLNOINC  MA(NET 
IN   x-I'LANE 

16.0 

4.0 

»/2  (ess) 

HALF- ATI  PTlfRE  OF 
HIDING  MAC  NET 
IN  y-PLANE   (jap) 

16.0 

S.O 

1/2  (cm) 

LENGTH  OF  SYSTEM 

16.0 

6.0 

L  (metres) 

FRINd   FIELD  0OR- 
RECTION  COEFFICIENT 

16.0 

7.0 

K|   (diAenstonless) 

FRINGE   F1F.I  D  TOR- 
RtCIION  UUcFFICIENT 

16.0 

1.0 

«,    (dtaenstonless) 

CURVATURE  OF 
ENTRANCE   FAlI  OF 
N.NUING  MAGNET 

16.0v 

12.0 

(1/R,)   (L/mctres) 

OPVATURE  OF 
EXIT   FACE  OF 
RLNUING  MAGNET 

16.0V 

13.0 

(1/R,)   (1/metres) 

FOCAL  PIANF. 
KHAIION 

16.0 

IS.O 

Angle  of  focal  plane  rotation   (degrees). 
See  type  code  16.0  for  details. 

INIT1AI.  BIAM  LINE 
x  -  COORD  1  NAP: 

16.0V 

16.0 

"• 

INITIAL  BEAM  LINE 

y-OD  RUINATE 

16.0V 

17.0 

r» 

INITIAL  BFAM  LINE 
Z-COORPINATE 

16.0V 

11.0 

«• 

INITIAL  BTAM  LINE 
lURI  KINTAL  ANGLE 

16.0V 

19.0 

•• 

INITIAL  BEAM  LINE 
VERTICAL  ANCLE 

16.0V 

20.0 

♦. 

SEOONn-ORTER 
CALCULATIONS 

17.0 

SEXTUPOLE 

IB.Ov 

LENCTH   (metres) 

FIELD  (kG) 

HALF-APERTURE  (ess) 

SOLENOID 

19. w 

LTNCTH    (metres) 

FIELO  (kC) 

BEAM  ROTATION 

20. v 

ANCLE  OF 

ROTATION  (degrees) 

STRAY   FIELD 

21.0 

See  later  section  of  report. 

Note:       The  v'l   foil 
Vary  Codes. 

owing  the  t 
The  units 

y-pe   codes   indicate   the  parameters  which  may  be  varied.      See   section  under   type  CO 
■te  standard  TRANSPORT  units   las  shown)  unless  changed  via  type  code   li.O  entries 

de  10.0  for  a  detalle 

d  explanation  of 

1 
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The  various  units  that  may  be  changed  are: 


Code 
Digit 


Quantity 


Standard 
TRANSPORT  Unit 


Symbols  used  in 
SLAC-75 


1.  0       horizontal  and  vertical  trans- 
verse dimensions,  magnet 
apertures,  and  misalignment 
displacements 

2.0       horizontal  and  vertical  angles 
and  misalignment  rotation 

angles 

*\ 
3.0       vertical  beam  extent  (only)  ' 

and  bending  magnet  gap  height 

*\ 

4.  0       vertical  beam  divergence   ' 

(only) 

5.  0       pulsed  beam  length  and  wave 

length  in  accelerator 

6.  0       momentum  spread 

7.0       bend,  pole  face  rotation,  and 
coordinate  layout  angles 

8.0       length  (longitudinal)  of  ele- 
ments, layout  coordinates, 
and  bending  magnet  pole  face 
curvatures 


cm 


mr 


cm 


mr 


cm 


percent  (PC) 
degrees  (DEG) 

metres  (M) 


x,y 


e,<P 


<P 


9.0 

magnetic  fields 

kG 

B 

10.0 

mass 

electron  mass 

m 

11.0 

momentum  and 

GeV/c 

P(0) 

energy  gain  in  ace 

elerator 

GeV 

AE 

section 

*)  These  codes  should  not  be  used  if  the  coordinate  rotation  (20.  0)  type  code  is 
used  anywhere  in  the  system. 
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